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Abstract 

Within local-spin-density functional theory, we have investigated the 'dissociation' of few- 
electron circular vertical semiconductor double quantum ring artificial molecules at zero magnetic 
field as a function of inter-ring distance. In a first step, the molecules are constituted by two 
identical quantum rings. When the rings are quantum mechanically strongly coupled, the 
electronic states are substantially delocalized, and the addition energy spectra of the artificial 
molecule resemble those of a single quantum ring in the few-electron limit. When the rings are 
quantum mechanically weakly coupled, the electronic states in the molecule are substantially 
localized in one ring or the other, although the rings can be electrostatically coupled. The effect 
of a slight mismatch introduced in the molecules from nominally identical quantum wells, or from 
changes in the inner radius of the constituent rings, induces localization by offsetting the energy 
levels in the quantum rings. This plays a crucial role in the appearance of the addition spectra as 
a function of coupling strength particularly in the weak coupling limit. 

PACS numbers: 85.35.Be, 73.21.-b, 73.22.-f, 71.15.Mb 
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I. INTRODUCTION 



Semiconductor quantum dots (QD's) are widely regarded as artificial atoms with prop- 
erties similar to those of 'natural' atoms. One of the most appealling is the capability of 
forming molecules. Systems composed of two QD's, quantum dot artificial molecules (QDM), 
coupled either laterally or vertically, have been investigated experimentally and theoretically 
at zero magnetic field (B), or submitted to magnetic fields applied in different directions, 
see e.g. Refs. MM , and references therein. 

Semiconductor ring structures have also received considerable attention in connection 
with the Aharonov-Bohm effect,— and the energy spectrum of nanoscopic self-assembled 
quantum rings occupied by few electrons has been experimentally analyzed. 7 Recently, high 
quality quantum rings (QR) have been fabricated on a AlGaAs-GaAs heterostructure con- 
taining a two-dimensional electron gas, by nano-litography with a scanning force microscope, 
see Refs. fa and |3, and references therein. These studies have allowed to extend previous 
studies to many-electron nanoscopic rings, and have provided an experimental determina- 
tion of the spin ground states of the rings by Coulomb-blockade spectroscopy, as well as the 
clear identification of a singlet-triplet transition, and the size of the exchange interaction 
matrix element,— properties that had been also determined in the past for QD's . 11 ^ 2 

Very recently, two different types of nanometer-sized QR complexes have been real- 
ized. One such complex consists of two concentric QR's grown by droplet epitaxy on an 
Alo.3Gao.7As substrate.— The other complex consists of stacked layers of InGaAs/GaAs 
QR's, whose optical and structural properties have been characterized by photolumines- 
cence spectroscopy and by atomic force microscopy, respectively.—^ Motivated by these re- 
cent experimental works, we have undertaken a theoretical study, within local- spin- density 
functional theory (LSDFT), of the ground state (gs) properties of QR's complexes at B = 0. 
In this work, we present the results we have obtained for the case of two GaAs vertical dou- 
ble quantum rings, leaving aside for a separate study the case of concentric double quantum 
rings, whose phenomenology is somewhat different.— To some extent, our work parallels the 
one we have carried out in the past for double QD's, with the aim of understanding the 
electronic properties of QRM's and the difference between vertical QDM and quantum ring 
molecule (QRM) structures. 

This work is organized as follows. In Sec. II we present the method we have used to 
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describe QR's and vertically coupled QRM's. The results we have obtained are discussed in 
Sec. Ill, and a brief summary is presented in Sec. IV. 



II. LSDFT DESCRIPTION OF QUANTUM RINGS AND VERTICAL QUANTUM 
RING MOLECULES 

We closely follow the method of Ref. Q, where the interested reader may find it de- 
scribed in some detail. We recall that within LSDFT, the ground state of the system is 
obtained by solving the Kohn-Sham equations. The problem is simplified by the imposed 
axial symmetry around the z axis, which allows one to write the single particle (sp) wave 
functions as (p n i a (r, z, 9, a) = u n i a {r, z)e~ tW Xa with n = 0, 1, 2, ...,/ = 0, ±1, ±2, . . ., being 
—I the projection of the sp orbital angular momentum on the symmetry axis. 

We have used effective atomic units h = e 2 /e = m =1, where e is the dielectric constant, 
and m the electron effective mass. In units of the bare electron mass m e one has m = m*m e . 
In this system, the length unit is the effective Bohr radius = aoe/m*, and the energy unit 
is the effective Hartree H* = Hm*/e 2 . In the numerical applications we have considered 
GaAs, for which we have taken e = 12.4, and m* = 0.067. This yields ~ 97.9 A and 
H* ~ 11.9 meV 

In cylindrical coordinates the KS equations read 



where r] a —+l(— 1) for cr=f(j), ^c/( r ? z ) i s the confining potential, V H (r, z) is the direct 
Coulomb potential, and V xc = d£ xc (n,m)/dn\ gs and W xc = d£ xc (n,m)/dm\ gs are the vari- 
ations of the exchange-correlation energy density S xc (n,m) in terms of the electron density 
n(r, z) and of the local spin magnetization m(r, z) = n\r, z) — n^(r, z) taken at the gs. 

As usual, £ xc (n,m) = S x (n,m) + S c (n,m) has been built from 3D homogeneous elec- 
tron gas calculations. This yields a well known, 18 simple analytical expression for the ex- 
change contribution £ x (n,m). For the correlation contribution £ c (n,m) we have used the 
parametrization proposed by Perdew and Zunger.— 




(1) 



+ V H + V xc + W xc r) a u n i a (r, z) = e nla u nla (r, z) , 
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For a double QR the confining potential V c f(r, z) has been taken parabolic in the xy plane 
with a repulsive core around the origin, plus a symmetric double quantum well of width w 
each, in the z direction. The potential in the xy plane has circular symmetry, and in terms 
of the cylindrical coordinate r it is written as 

V cf (r) = Vo&(R -r) + ^muj 2 (r-R ) 2 e(r-R ) , (2) 

with 0(x) = 1 if x > and zero otherwise. The convenience of using a hard- wall confining 
potential to describe the effect of the inner core in QR's is endorsed by several works in the 
literature.-^ We have taken Rq = 5 nm, w = 5 nm, Vq = 350 meV, and u = 15 meV. The 
depth of the double quantum well is also Vq. This set of parameters fairly represents the 



smallest rings synthesized in Ref. |22j, and together with the distance d between constituent 
quantum wells, determine the confining potential. The distance d is varied to describe 
quantum ring molecules at different inter-ring distances. For the single 'thick' QR we will 
discuss as a reference system, we have used the same confining potential in the xy plane, 
together with a single quantum well in the z direction. For all structures, the sharp potential 
wells have been slightly rounded off, as shown in Ref. |3- Details about how the KS and 
Poisson eqs. have been solved can be found in Ref. [l?! 



III. RESULTS AND DISCUSSION 



A. Single quantum ring 

We have carried out calculations for a single thick QR confined as indicated in the previous 
Section, and for a strictl y tw o-dimensional QR confined by the radial potential Eq. (J2J), as 



indicated, e.g., in Ref. |23|. These results will help us to discuss the appearence of the 
addition spectra of the QRM. 

Fig. ^ shows the addition energies A 2 (N) 



A 2 (N) = E(N + 1) - 2E(N) + E(N - 1) , (3) 

where E(N) is the total energy of the N electron QR, as a function of N. It can be seen that 
the 2D and 'thick' -i.e., axially symmetric 3D- models sensibly yield the same results for 
this observable, a well know result for QD's.i For the thick rings, the value of the calculated 
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total spin third component, 2S Z , is also indicated in the figure. We want to point out that 
in the N = 3 case, the 2D model configuration is fully polarized (2S Z = 3). This is due 
to the fact that the exchange-correlation energy is overestimated by strictly 2D models.— 
Fully polarized N = 3 QR configurations are not an artifact of the LSDFT. As a matter of 
fact, they have been also found by exact diagonalization methods for some ring sizes and 
confining potential choices. 25 



The gs spin assignments we have found here coincide with those of Ref. |23j, although the 
height of the peaks in A 2 (iV) depends to a large extent on the confining potential. They 
are related to the relative stability of the electronic shell closures in the ring, which for 
N > 6 are substantially different from these of QD's. In the case of rings, they are mainly 
governed by the fourfold degeneracy of the non- interacting sp levels with \l\ ^ 0, and the 
twofold degeneracy of the non-interacting sp levels with |/| = 0. This yields the marked shell 
closures at iV = 2, 6, 10, 20 and 28 with S z = 0, as well as the S z = gs found for N = 24. 
The 2S Z = 2 ground states that regularly appear between them indicate that Hund's rule is 
fulfilled by single QR's. 

The complex spin structure around N = 13 deserves some comments. It is due to the 
occupancy of the second (s |) state with I = -this spin structure is missing in other QR 
calculations that employ a different confining potential.— Figure |2 displays the sp energies 
e n \ a for N = 13, which are distributed parabolic-like as a function of I, each parabola 
corresponding to a different value of the principal quantum number n. This figure explains 
the filling sequence around N = 13. For N = 12, the second (0 |) state is empty, yielding 
2S Z = 2; for N = 13 the exchange interaction favors the filling of this state yielding 2S Z = 3; 
for iV = 14, one of the (±3 j) states is filled -they are degenerate-, yielding 2S Z = 2 (actually, 
this many-electron configuration is nearly degenerate with the one in which the (0 J.) state 
is filled instead, which also yields 2S Z = 2). For N = 16, the (0 j) and (±3 j) become 
populated, producing a fairly strong shell closure. 



B. Homonuclear quantum ring molecules 

We consider first the case of a QRM formed by two identical quantum rings. By anal- 
ogy with natural molecules, we call them homonuclear QRM. We have calculated their gs 
structure for d — 2, 4 and 6 nm, and up to N = 32. For a given electron number N, the 
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evolution of the gs ('phase') of a QRM as a function of d may be thought of as a dissocia- 
tion process.— Within LSDFT, each sp molecular orbital has, as quantum labels, the third 
component of the spin and of the orbital angular momentum, the parity, and the value of 
reflection symmetry about the z = plane. Symmetric states \S) are called bonding states, 
and antisymmetric states \AS) are called antibonding states. 

The energy splitting between bonding and antibonding sets of sp states, Asas, can be 
properly estimated^ from the energy difference of the antisymmetric and symmetric states 
of a single electron QRM, A S as ~ -^( 2 ^« ) — -^( 2 ^<j~ ) ~ see below for the notation-, and varies 
from 24.9 meV at d = 2 nm (strong coupling), to 1.49 meV at d = 6 nm (weak coupling). 
In this range of inter-ring distances, As as can be fitted as As^s = A e~ d ^ d °, with A = 82 
meV and do = 1.68 nm. The relative value of the two energies Kujq and A5/L5 crucially 
determines the structure of the molecular phases along the dissociation path. 

Figure H3 shows the evolution with d of the gs energy and molecular phase of a QRM 
made of N = 3 — 7 electrons. Each configuration is labeled using an adapted version of the 
ordinary spectroscopy notation,— namely +1 L* where S is the total \S Z \, and L is the 
total \L Z \. The superscript +(— ) refers to even (odd) states under reflection with respect to 
the z = plane, and the subscript g(u) refers to positive (negative) parity states. To label 
the molecular sp states we have used the standard convention of molecular physics, using 
a, 7r, 6, . . . if I = 0, ±1, ±2, . . .. Upper case Greek letters are used for the total \L Z \. Fig. |3] 
shows that the energy of the molecular phase increases with d. This is due to the increase 
of the energy of the sp bonding states as d increases,-^ that dominates over the decrease 
in Coulomb energy. At larger inter-ring distances, the constituent QR are so apart that 
eventually the decrease of Coulomb energy dominates and the tendency is reversed. The 
phase sequences are the same as for double quantum dots,— although the transition inter-ring 
distances, which obviously depend on the kind and strength of the confining potential, are 
different. As for double quantum dots, we have found that the first phase transition of a 
few-electron QRM is always due to the replacement of an occupied bonding sp state by an 
empty antibonding one. 

Figure 0] shows the addition spectra for homonuclear QRM up to N = 31 for the three 
selected inter-ring distances. Also shown is the reference spectrum of a single QR. The 
spectra have been offset for clarity. For small d (As as ^ ^0) the spectrum of the QRM 
is rather similar to a single QR, especially for few-electron systems, with minor changes 
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arising in the N ~ 12 and ~ 24 regions that will be commented below. It is clear that for 
d = 2 nm the two QR's are electrostatically and quantum-mechanically coupled and behave 
as a single system. At intermediate distances the spectrum pattern becomes more complex, 
but at larger distances (e.g. d = 6 nm), when the QRM molecule is about to dissociate, 
the physical picture that emerges is rather simple and can be interpreted using intuitive, 
yet approximate arguments: At large distances (As^s <C hu ), the QR's are coupled only 
electrostatically, and most \S) and \AS) states are quasidegenerate. Electron localization in 
each constituent QR can be achieved combining these states as (\S) ± \AS))/y2 and as a 
consequence, the strong S z = peaks found at N = 12 and 20 are readily interpreted from 
the peaks appearing in the single QR spectrum at iV = 6 and 10; the process can be viewed 
as the symmetric dissociation of the original QRM leading to very robust closed shell single 
QR configurations. This is also the origin of the QRM S z = peaks at iV = 2 and 4. In the 
former case, the QRM configuration corresponds to one single electron being hosted in each 
constituent QR coupled into a singlet state, and in the latter case, the QRM configuration 
is viewed as two QR's, each one occupied by two electrons filling the Is shell. 

At d = 6 nm, other dissociations display a more complicated pattern, such as 16 — > 8 + 8, 
or 8 — > 4 + 4, whose final products are QR's that fulfill Hund's rule whereas the actual QRM 
has S z = 0. These could be interpreted as rather entangled QRM, 'harder' to dissociate, for 
which a d — 6 nm inter-ring distance is not large enough to allow for electron localization. 
The fact is that not only quasidegeneracy of occupied \S) and \AS) states at given d plays a 
role in this intuitive analysis, but also whether their number is equal or not, so that they may 
be eventually combined to favor localization. An example of these two different situations 
is illustrated in Fig. where we show the sp states of the N = 16, 20 and 23 QRM at 
d = 6 nm. In the case of N = 16 and 23, the filled bonding states near the Fermi level 
have not filled antibonding partner and are delocalized in the whole volume of the QRM, 
contributing to the molecular bonding at that distance, whereas all other bonding states 
can be localized combining them with their antibonding partner: as in natural molecules, 
some orbitals contribute to the molecular bonding, whereas some others do not. 
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C. Heteronuclear quantum ring molecules 

For vertically coupled lithographic double quantum dots, it has been found unavoidable 
that a slight mismatch is unintentionally introduced in the course of their fabrication from 
materials with nominally identical constituents quantum wells J« which is responsible for 
electron localization as the interdot coupling becomes weaker. This offsets the energy levels 
in the quantum dots by a certain amount that was there estimated to be up to 2 meV, and 
this plays a crucial role in the appearance of the addition energy spectra as a function of the 
coupling strength particularly in the weak coupling limit. A similar picture is also found in 
coupled self-assembled quantum dots, where strain propagation between adjacent layers of 
dots often leads to top QD's of increased size.— 

Likely, the same fabrication limitations will appear in the case of vertically coupled dou- 
ble quantum rings. Anticipating to this situation, we have carried out a series of QRM 
calculations in which the double quantum wells have the same width w but slightly different 
depths, namely Vb±5, with 5 <C Vq. It can be easily checked that in the weak coupling limit 
(25 3> &sas)i 25 is approximately the energy splitting between the bonding and antibonding 
sp states, which would be almost degenerate if 5 = 0. For this reason we call the mismatch 
(offset) the quantity 28. 

We have considered two possible values of the mismatch, namely 25 = 2 and 4 meV, and 
have obtained the corresponding addition spectra for up to N = 13 electrons -according to 
our previous experience with double quantum dots,- we expect that the larger differences 
will arise in few-electron QRM. The results are displayed in Fig. |U1 It can be seen that in 
the strong coupling limit, the effect of the mismatch on the addition energies is negligible, 
as expected.— The electrons are completely delocalized in the whole volume of the QRM, 
and the introduced mismatch is unable to localize them in either of the constituents QR's. 

For the few-electron QRM, which is the more interesting physical situation, we have 
shown before that the fingerprint of homo nuclear character is the appearance, in the weak 
coupling limit, of the peaks in the addition spectrum corresponding to N = 2 and 4, as well 
as the spin assignment S z = 0. It can be seen from Fig. El that in the intermediate regime 
(d = 4 nm) the N = 4 peak still corresponds to a 2S Z = 2 configuration, but at larger inter- 
ring distances, it eventually disappears, yielding an addition spectrum that clearly manifests 
the heteronuclear character of the QRM and constitutes a clean fingerprint of these kind of 
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configurations. 

It is useful to display the dissociation of the QRM representing the d evolution of the sp 
molecular wave functions, introducing the ^-probability distribution function 1 



Examples of these probability functions can be seen in Fig. [7[ where we show V(z) for 
iV = 20, 25 = 4 meV, and iV = 8, 25 = 2 meV (deeper well always in the z < region), 
each for the chosen three d values. In each panel the probability functions are plotted 
ordered from bottom to top according to the increasing sp energies. For N = 20 the final 
configurations are the closed shell iV = 10, 2S Z = QR's, whereas for N = 8 the A = 4, 
2S Z = 2 Hund's rule QR configurations emerge. 

Finally, we discuss the case of two QR's of different radii vertically coupled to build an 
axially symmetric QRM, and study the effect this assymetry has on the addition spectrum 
(we have discarded a possible disalignment of the QR symmetry addressing this 

situation would require a much more demanding full 3D calculation 28 ). To this end, we have 
taken for one ring R = 6 nm, while for the other one we have kept the same value as before, 
Ro = 5 nm (5 is set to zero this case). 

We show in Fig. [HJ the addition spectra for up to N = 14 electrons and d — 2, 4, and 6 
nm. It can be seen that in the strong and intermediate coupling cases they are fairly similar 
to the previous heteronuclear case -and to the homonuclear case as well-, indicating a fairly 
robust structure of the QRM in these limits. As before, the heteronuclear character clearly 
shows up in the weak coupling limit, with a peak structure and S z assignments remarkably 
similar to those discussed in the previous situation with 5^0. 

IV. SUMMARY AND OUTLOOK 

We have discussed the appearance of the addition energy spectra of homonuclear and het- 
eronuclear quantum ring molecules at zero magnetic field. In particular, we have addressed 
the addition energy spectrum of QRM from the weak to the strong coupling limits. Finger- 
prints of homo- and heteronuclear molecular character have been pointed out in the weak 
coupling limit. As it happened in the study of vertically coupled double quantum dots,— we 
believe this may be helpful in the analysis of future experiments on vertical QRM's. 




(4) 
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The present study can be naturally extended to the case of QRM's submitted to magnetic 
fields of arbitrary direction. A rich interplay between molecular phases having different 
isospin is expected to appear as a function of which might have an observable influence 
on the Aharonov-Bohm effect and on the far-infrared spectroscopy of nanoscopic QRM's. 
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FIG. 1: Addition energies A2(N) (meV) as a function of the number of electrons N for a thick 
QR (solid dots, solid lines) and a strictly 2D QR (open dots, dashed lines). The value of 2S Z is 
indicated. 
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FIG. 3: Energy (meV) and gs molecular phases of the homonuclear QRM as functions of the 
inter-ring distance d for N = 3 — 7. For each QRM, different phases are represented by different 
symbols. 
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FIG. 4: A2^N) for homonuclear QRM with inter-ring distances d = 2,4, and 6 nm. The addition 
energies have been offset for clarity. Also shown is the reference spectrum for a single QR. The 
value of 2S Z is indicated. 
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FIG. 5: Single particle energy levels (meV) as a function of I for an homonuclear QRM with 
d = 6 nm and iV = 16 (top panel), N = 20 (middle panel), and iV = 23 (bottom panel) . 
Upward (downward) triangles denote f(J,) spin states. Open(solid) triangles correspond to anti- 
bonding(bonding) states. The thin horizontal line represents the Fermi energy. 



15 




FIG. 6: A2(-/V) for heteronuclear QRM with inter-ring distances d = 2,4, and 6 nm and mismatch 
25 = 2 meV (left panels) and 25 = 4 meV (right panels). The addition energies have been offset 
for clarity. Also shown is the reference spectrum for a single QR. The value of 25*2 is indicated. 
Note that in some cases two different values of S z have been asigned to the same peak. This means 
that the corresponding configurations are nearly degenerated. 



16 



N=20 25=4meV 



d=2nm £ 


d=4nm £ 


d=6nm L 


±jiU 








. I A 


atl 
/V\ 




**Aa 




A ±814 


±4>TI 
/V\ 




a!^ 




±nti 


±8ti 
/V\ 




aA 






±jiU 
/V\ 




A^ 






atl 

. /V\ . 








A A, 



-15 



15 -15 



15 -15 15 



N=8 28=2meV 



d=2nm £ 



±8t 
M 

±7lT 



d=4nm £ 



d=6nm E 



A A 
±7iT4- 



oU 
A/\. 



±Jti 
/\ 


A 


A 


±7tt 




A 


/\ 


o-l 


A 


A 





-15 15 -15 15 

z (nm) 



■15 15 



FIG. 7: Calculated probability distributions V[z) (arbitrary units) as a function of z for heteronu- 
clear QRM's with TV = 20, 25 = 4 meV (top panels), and TV = 8, 25 = 2 meV (bottom panels). 
The corresponding molecular configuration is also indicated. 
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FIG. 8: A2(iV) for heteronuclear QRM made of QR with different core radii Rq = 5 and Rq = 6 
nm, hosting up to N = 14 electrons and inter-ring distances d = 2,4, and 6 nm. The addition 
energies have been offset for clarity. Also shown is the reference spectrum for a single QR. The 
value of 2S Z is indicated. 
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